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Abstract 

The assumption of collisionless cold dark matter on its own cannot reconcile several astro- 
physical discrepancies (cusp-vs-core problem, missing satellite problem, too-big-to-fail problem). 
Self-interacting dark matter provides a promising framework for solving all these problems, and 
self-interaction cross sections are duly constrained in the literature. Following the work of Tulin, 
Yu, and Zurek [5, we can constrain the dark matter mass and the mass of a light mediator as¬ 
suming a generic scalar Yukawa-type interaction. In particular, we constrain the strongly coupled 
inflationary dark matter of the luminogenesis model, a unification model with the gauge group 
SU{3)c X SU{6) X U{1)y, which breaks to the Standard Model with an extra gauge group for 
dark matter when the inflaton rolls into its true vacuum. The luminogenesis model is addition¬ 
ally subject to constraints on inflation, and we find an upper bound on the scale of symmetry 
breaking of the inflaton and the decoupling scale Mi of certain representations of the gauge group. 
We emphasize that the luminogenesis model enables a unique connection between astrophysical 
constraints, the nature of dark matter, and inflation. 
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Introduction 


The formation of galaxies and galaxy clusters is heavily influenced by the nature of dark 
matter. For the usual framework of cold dark matter, there are discrepancies between their 
predictions and observations. A^-body simulations for exclusive collisionless cold dark matter 
predict the central density prohle of dwarf galaxy and galaxy cluster halos to be very cusp- 
like, whereas observations indicate flat cores (cusp-vs-core problem) [2]. The number of 
Milky Way satellites predicted in simulations is bigger by an order of magnitude than the 
number inferred from observations (missing satellite problem) [3l 0], although this may not 
be very troublesome if more ultra-faint galaxies are successfully detected in the future [5]. 
The brightest observed dwarf spheroidal galaxy satellites of the Milky Way are predicted to 
be in the largest Milky Way subhalos, but the largest subhalos are too massive to host them 
(too-big-to-fail problem) [6]. The resolution of these problems may come through several 
possible means, including more accurate consideration of baryon interactions, astrophysical 
uncertainties, and warm dark matter. A promising framework that can solve all these issues 
is self-interacting dark matter. 

In this work, we are interested in astrophysical constraints on strongly coupled scalar 
Yukawa-type interactions, and our methodology follows that of pQ. (And we explain our 
use of the term ”Yukawa-type” later in the text.) The astrophysical constraints on self¬ 
interaction cross sections confer restrictions on the allowed masses of dark matter and the 
scalar light mediator in the Yukawa-type interaction. When applied to inflationary dark 
matter of the luminogenesis model m El □□], these restrictions imply constraints on the 
unihcation scale of the model and the parameters of the inflation potential. 

We first briefly review the luminogenesis model and its significance. We then present 
our results constraining the dark matter and scalar mediator masses that come from com¬ 
paring the Yukawa-type interaction cross section from solving Schrodinger’s equation and 
astrophysical constraints. Applying the results to the luminogenesis model, we arrive at 
a unique, interesting connection between astrophysical bodies, dark matter, and cosmic 
inflation. 


Review of the Luminogenesis Model 

In the luminogenesis model, the dark and luminous sectors are unihed above the Dark 
Unihed Theory (DUT) scale. At this DUT scale, the unified symmetry of the model breaks 
{SU{3)c X SU{Q) X U{1)y ^ SU{3)c x SU{4:)dm x SU{2)l x U{1)y x U{1)dm), and the 
breaking is triggered by the inflaton’s slipping into the minimum of its symmetry-breaking 
(Coleman-Weinberg) potential and acquiring the true vacuum expectation value iidut, which 
is the DUT scale energy. This symmetry breaking allows the inflaton to decay to dark 
matter, and dark matter can in turn decay to Standard Model (SM) and ’’mirror” matter. 
The representations and group structure of the luminogenesis model for each of the three 
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families are given below [8]. The existence of ’’mirror” fermions, as proposed by 0 , is 
necessary for anomaly cancellation, and it provides a mechanism in which right-handed 
neutrinos may obtain Majorana masses proportional to the electroweak scale, and they 
could be searched for at the Large Hadron Collider. 


SU{Q) 

SU{A)dm X SU{2)l X U{1)dm 

6 

(1,2; 

2 +(4,l)-i 

20 

(4,1' 

3 + (4L1)_3 + (6,2)o 

35 

(1,1' 

0 + (15,1)0 + (1, 3)0 + (4, 2)_3 


* 

,2)3 


TABLE I. ( 1 , 2)2 represents luminous matter while ( 4 ,1)3 -|- ( 4 *, l)-3 represent dark matter. 



5C(3)c X SU{6) X U{1)y 

R D SM fermions 

(3,6,1/6)l + (1,6,-1/2)l 
+(3,1,2/3)r + (3,1,-1/3)r 
+(1,1,-1)r 

R D ’’mirror fermions” 

(3,6,1/6)r + (1,6,-1/2)r 
+(3,1,2/3)l + (3,1,-1/3)l 
+(1,1,-1)l 

R D dark-matter fermions 

(1,20,0) 


TABLE II. R in the left column denotes representation. Standard Model (SM) left-handed doublets 
and right-handed singlets comprise the first row, mirror right-handed doublets [ 9 l [TO] and left- 
handed singlets comprise the second row, and dark-matter left- and right-handed fermions belong 
to the last row. 

Mirror quarks and leptons can be searched for at the LHC m- The search for the 
electroweak-scale right-handed neutrinos is particularly interesting since it will be a direct 
test of the seesaw mechanism. As emphasized in [S] , the production of i^RiyR at the LHC can 
give rise to interesting signals such as like-sign dileptons. Furthermore, the lightest mirror 
quark can decay into a Standard Model (SM) quark by emitting a SM-singlet Higgs scalar 
(Qr Ql + (j^s) through an interaction Lagrangian of the form gsqQL (j^s Qr + h.c., where 
ql and refer to a SM left-handed and mirror right-handed quark doublet respectively. A 
similar decay process applies to the lightest mirror lepton (/^ —)■ Il + 4>s) with gsi Ir 4>s ^r + 
h.c. 

The inflaton is represented by (1, l)o of 35, and since 20 x 20 = Ig -l- 35a + 17’5s -l- 189a, 
the inflaton decays mainly into dark matter through the interaction g 2 Q 'Ljo'^ 2^20 </’35 which 
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contains the inflaton in g2oX^^2X^4>inf- This is one of the main points of [ 8 ]: the predomi¬ 
nance of dark over luminous matter. Dark matter can then be converted via the process of 
what [H] refers to as ’’luminogenesis.” We summarize this process below: 

• As is noted in [S], for the SU{‘i)]:,M dark matter (DM) fermion x, a small asymmetry 

in the number density = n^—n^ is assumed to be present, with +An^ 

only slightly bigger than = ngym, where Usym is the symmetric part of and 
n^. The asymmetric part consists of the small excess An^ Usym- For the origin of 
this asymmetry in the DM number density, we assume that there is a global U{1)^ 
symmetry for DM. The interactions involving the gauge bosons of the coset group 
SU ( 6 )/ SU (4) X SU (2) X f/(l) dm explicitly break the ^(1)^ symmetry, and their decays 
involving the interferences between the tree-level and one-loop diagrams will ultimately 
produce a net DM asymmetry, assuming the presence of CP violation in the DM sector. 
(This process is similar to the one involving X and Y gauge bosons in SU{5) Grand 
Unihed Theory.) 

• X and X can annihilate via xdm, the massive dark photon of U{1)dm, into luminous 
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particle-antiparticle pairs via the effective interaction — {xiu.x){fl^f)i nnd the 
particle-antiparticle pairs of luminous fermions annihilate to radiation. 

• X nnd X can also be converted to luminous leptons via the interactions with two 

scalar helds: $^ 5 ^ and and given by {Xl'^ 21 l) + h.c., resulting in 

Xl + Xr Il+Ir and Xl+Xr ^l + Ir ■ Although the decay length of the previously- 
mentioned process Ir ^ II + (j^s could be macroscopic at the LHC (’’long-lived” l^), 
in the early universe, basically decays promptly into SM leptons. 

• The coefficients for the annihilation process involving Xdm and the conversion process 
involving the scalars $[ 5 ^ and are independent of each other, and they are such 
that 14% of all dark matter (14% of asymmetric and symmetric parts) converts via the 
scalars and 86 % of asymmetric and symmetric parts annihilates via Xdm ultimately to 
radiation. The 14% of the symmetric parts of and that is converted to luminous 
matter via the scalars has equal parts of luminous particles and anti-particles, and 
these annihilate to radiation. So the whole of the symmetric parts of and is 
mainly converted to radiation. But since annihilation via Jdm requires the presence of 
both X and x, this annihilation does not affect the asymmetric part An^. So overall, we 
are left with radiation from the symmetric parts of and n^, and of the asymmetric 
part, 14% is luminous matter and 86 % is DM, giving the correct proportion of luminous 
to dark matter as the conversion process via the scalars freezes out. 

So we see that this process of luminogenesis depends upon both freeze-out and asymme¬ 
try in dark matter (and the asymmetry is propagated through to luminous matter), and the 


3 



coefficients of the processes discussed are such that luniinogenesis gives what is observation- 
ally expected. These processes involved in luniinogenesis are discussed in more detail in [ 8 ] 
and are being explored further in |12j . 

It is assumed that (15, l)o + (l,3)o -|-(4, 2 )_ 3 -|- (4*, 2)3 of 35 and (6, 2)o of 20 have 
masses that are on the order of the DUT scale and thus do not affect the particle theory below 
that energy scale. The SU{4:)dm dark matter fermions are represented by (4, 1)3 -|- (4*, 1)_3 
in the 20 representation of SU{6). Since dark matter should have no U{1)y charge, the 
SU{4:)dm particles in (4, l)_i in the 6 representation of SU{6) cannot be dark matter since 
they have t/(l)y charge (as shown in Table 0 . and they are assumed to decouple below the 
mass scale we call Mi. In [13], we make predictions for the mass of y in the following way: 

• We run the SU{2)l gauge coupling from the known electroweak scale up to some 

unknown DUT scale where it intersects with the gauge coupling a^. 

• Then we run 0:4 down to its conhnement scale, which is when 0:4 ~ 1. In analogy with 
QCD conhnement of SU{3)c, the main contribution to S't/(4)/)M fermions’ dynamical 
mass is from the conhnement scale of S'f/(4)£)M, and that energy scale is our mass 
prediction for y. 

• In order to specify that scale, we need to specify a DUT scale. Since SU{6) breaks 
at the DUT scale when the inhaton slips into its true vacuum, we specify the DUT 
scale and therefore the dynamical mass of y by constraining the parameters of a 
symmetry-breaking (Coleman-Weinberg) inhaton potential with Planck’s constraints 
on the scalar spectral index and amplitude nans]- 

We will compare these predictions from [T3] to the results from this paper’s application of 
astrophysical constraints to dark matter mass. 

Because of the conhnement of S'f/(4), dark baryons are formed from four y particles. 
These particles are dubbed CHIMPs, which stands for ”y Massive Particles.” A CHIMP 
is denoted by X, and X = (yyyy). As we know from QCD, SU{3) Nambu-Goldstone 
(NG) bosons appearing from the spontaneous breaking of chiral symmetry from < qq 0 
acquire a small mass from the explicit breaking of quark chiral symmetry due to the small 
masses of quarks, and they become pseudo-NG bosons known as pions. We expect a similar 
phenomenon from < yy > 7 ^ 0 in 517(4), and the NG bosons can acquire a small mass 
through a term moyy with mo a parameter that is related to and should be on the 

order of and mo -C A 4 ~ We seek to constrain the m^^^j^-mx parameter space 

through astrophysical constraints via the procedure in the following section. 

Solving Schrodinger’s Equation 

For unspecihed X and ttdm, ffi general, the cross section of their interaction may not 
he in the regimes of the Born or classical approximations, so we cannot rely solely on 
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analytical expressions for these regimes. In order to find how the mass of strongly conpled 
DM is correlated to the mass of a scalar mediator via astrophysical constraints, we need to 
nnmerically solve Schrodinger’s eqnation, and we nse the methodology described in detail 

in d]. 

We take the interaction between DM (a CHIMP, denoted by X = (xXXX)) a scalar 
mediator {tidm) to be given by an attractive Yukawa-type potential 

V-(^) = (1) 

r 

with interaction 

= QOM^Xt: DM-, ( 2 ) 

where aoM = 9dm /{^^)- We use the term ’’Yukawa-type” because a Yukawa interaction 
is typically thought of as between a spin-t particle and a scalar, whereas we are using Eq. 
Q to describe the interaction between X, a boson, and the scalar tidm- This distinction 
does not make a difference in calculation in our case since Schrodinger’s equation does not 
account for spin. In analogy with QCD’s strong nucleon-pion interaction, we assume that 

<^DM ~ 1- 

For carrying out the computational method for solving Schrodinger’s equation exactly 
as described in [1] with a similar level of computational precision, we make plots of mx 
vs for aoM = 1 and udm = ‘^ via their relationship through the velocity-averaged 

transfer cross section < ar > for the interaction described by the potential in Eq. Q. The 
plots are given in Figs. Q and (|^. 

Using the convention of [1], the plots are described as follows: 

• Green lines going from left to right respectively represent {(Jt)/^x = 10 and 0.1 on 
dwarf scales, required for solving small scale structure anomalies. 

• Red lines going from left to right respectively represent {(Tt)I mx = 1 and 0.1 on Milky 
Way (MW) scales. 

• Green lines going from left to right respectively represent (crr)/'mx = 1 and 0.1 on 
cluster scales. 

The above astrophysical upper and lower bounds on (crT)/mx, as discussed in [1], come 
largely from N-body simulations for a limited number of specihc cross sections, so their 
constraining power in our plot should not be taken to be extremely stringent. But the 
ranges given for (cTr)/mx are generally what is needed to satisfy observational constraints 
from structure formation, and we discuss the regions of mx-m-^oM parameter space that fall 
within all three ranges (within the bounds of all three sets of colored lines) of {aT)/mx- 


Analysis of Results 
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For aoM ~ 1, we can see from Figs. Q and (|^ that all constraints from cluster, the 
Milky Way, and dwarf galaxies can be met for mx ranging from a few 100 GeV to a few 
TeV, and this range corresponds to 1 MeV < < 10 MeV. For almost all of mx < 100 

GeV, we found regions of parameter space that did not grossly disagree with the three 
astrophysical constraints, similar to the results of the resonant regimes in Fig. 6 of |1]. We 
focus our plots on the noteworthy observation that mx ^ 10 TeV does not agree with all 
three constraints in the plots (barring the fact that the tightness of these constraints is open 
to interpretation). 

Given the numerical results in the previous paragraph, and since A 4 ~ m^ < mx/4, one 
can see that the approximation mo ~ ^ ^4 seems to be a good one, much better than 

the analogous chiral approximation in QGD. This connection between the constraints on the 
macroscopic astrophysical scale and the microscopic ttom-X interaction lends support to the 
viability of the luminogenesis model. 



FIG. 1. We plot mx vs for the case of aoM = 1- We see that all three constraints from 

clusters, the Milky Way, and dwarf galaxies (described in the text) can be met for mx ranging 
from a few 100 GeV to a few TeV since this parameter space falls within all three sets of colored 
lines. 

We now consider the implications of this upper bound on the mass of strongly coupled 
DM for the luminogenesis model. Since we saw that X = (xXXX) cannot have a mass bigger 
than a few TeV, and since m^ < mx/4, we see there is an upper bound of 1 TeV for m^. 
In Fig. (g, we plot jJ-DUT vs Ml for this constraint m^ < 1 TeV using Eq. 10 derived from 
the unihcation of SU{2)l and SU{4:)dm at the DUT scale in PI- From Fig. (|^, we see 
that ^dut < 10^® GeV in order for this astrophysical upper bound for m^ to be satisfied, 
and most of the viable parameter space (the shaded triangle) is for values of Hdut much 
less than 10^® GeV. Along with this constraint, we also see that Mi < 10® GeV. 
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FIG. 2. We plot mx vs for the case of aoM = 2. We see that all three constraints from 

clusters, the Milky Way, and dwarf galaxies (described in the text) can be met for mx ranging 
from a few 100 GeV to a few TeV since this parameter space falls within all three sets of colored 
lines. 



Using this upper bound on hdut along with Planck’s constraints on the scalar spectral 
index and amplitude, we can also determine upper bounds on the number of e-folds and the 
parameters of the potential for inflation (in our case, the Coleman-Weinberg potential we 
used in |I3]). We work out the relationships of these parameters under the constraints from 
Planck in Eq. 21 of [13], and one can see that the number of e-folds would need to be less 
than roughly 95. 


7 














Conclusion 


In this work, we apply astrophysical constraints from clusters, dwarf galaxies, and the 
Milky Way obtained from N-body simulations for self-interacting dark matter in order to 
solve the structure formation problems known as the cusp-vs-core problem, the missing 
satellite problem, and the too-big-to-fail problem. Using the methodology of [1], we assume 
a Yukawa-type interaction between dark matter and a light mediator for a strongly coupled 
interaction. In particular, we examine the luminogenesis model and the interaction between 
the bosonic and baryonic CHIMP of S'17(4), X = (xXXX)^ and the light mediator tidm from 
chiral symmetry breaking. 

From the astrophysical constraints, for mx > 100 GeV, we saw that mx can range 
from a few 100 GeV to a few TeV, and this range corresponds to 1 MeV < fn-woM ~ 10 
MeV. These values imply a much better chiral approximation tuq -C A 4 than the analogous 
approximation in QGD. From our analysis, the upper bound of mx < a few TeV led to 
the upper bounds ^dut < 10^® GeV and Mi < 10® GeV. Using the upper bound on hdut 
along with Planck’s constraints on inflation, we can obtain upper bounds for the parameters 
of an inflation potential and the number of e-folds. For our Goleman-Weinberg model, the 
number of e-folds would need to be less than roughly 95. 

With the application of structure formation constraints on the self-interacting dark matter 
of the luminogenesis model, we have found an interesting observational connection between 
the macroscopic and the microscopic. We have established a unique connection between 
observational astrophysics, the nature of dark matter, and inflation, and we hope to exploit 
more connections in the future in order to further sharpen our understanding of dark matter 
and the origins of our universe. 


Acknowledgements 

PQH was supported in part by the U.S. Department of Energy under Grant No. DE- 
FG02-97ER41027 and is supported by the Pirrung Foundation. KJL is supported by the 
Pirrung Foundation. 


[1] Sean Tulin, Hai-Bo Yu, and Kathryn M. Zurek, Phys. Rev. D 87, 115007 (2013). 
arXiv:1302.3898[hep-ph]. 

[2] W. G. J. de Blok, Advances in Astronomy, 2010, 789293 (2010). 

arXiv:0910.3538[astro-ph.CD]. 

[3] Ben Moore et al, Astrophys. J. 524, L19 (1999). arXiv: astro-ph/9907411, 

[4] Anatoly A. Klypin et al, Astrophys. J. 522, 82 (1999). arXiv:astro-ph/9901240. 

[5] Joshua D. Simon and Marla Geha, Astrophys. J. 670, 313 (2007). 

arXiv:0706.0516[astro-ph]. 

[6] Michael Boylan-Kolchin, James S. Bullock, and Manoj Kaplinghat, Mon. Not. Roy. Astron. 
Soc. 415 ,L40 (2011). 1103.0007 [astro-ph.CO]. 

[7] P.H. Frampton and P.Q. Hung, Phys. Lett. B675, 411 (2009). arXiv:0903.0358[hep-ph]. 

[8] P.H. Frampton and P.Q. Hung, arXiv: 1309.1723 [hep-ph]. 

[9] P.Q. Hung, Phys. Lett. B649, 275, (2007). arXiv:hep-ph/0612004, 

[10] Vinh Hoang, P.Q., Hung, and Ajinkya Kamat, Nucl. Phys. B 10, 002 (2013). 

arXiv:1303.0428[hep-ph]. 

[11] S. Chakdar, K. Ghosh, V. Hoang, P. Q. Hung, T. Le, and S. Nandi, in preparation. 

[12] P. Q. Hung, work in progress. 

[13] P.Q. Hung and Kevin J. Ludwick, JGAP 09, 031 (2015). arXiv: 1411.1731 [hep-ph]. 

[14] P.A.R. Ade et al. (Planck Collaboration), Astron.Astrophys. 571, A16 (2014). 
arXiv:1303.5076[astro-ph.CD]. 

[15] P.A.R. Ade et al. (Planck Collaboration), Astron. Astrophys. 571, A22 (2014). 
arXiv:1303.5082[astro-ph.CD]. 


9 




